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ABSTRACT 

We present new measurements of the Ha luminosity function (LF) and star formation rate (SFR) 
volume density for galaxies at z ~ 0.8. Our analysis is based on 1.18/im narrowband data from 
the NEWFIRM Ha (NcwHa) Survey, a comprehensive program designed to capture deep samples 
of intermediate redshift emission-line galaxies using narrowband imaging in the near-infrared. The 
combination of depth (w 1.9 x 10~ 17 erg s" 1 cm" 2 in Ha at 3tr) and areal coverage (0.82 deg 2 ) of 
the 1.18/im observations complements other recent Ha studies at similar redshifts, and enables us to 
minimize the impact of cosmic variance and place robust constraints on the shape of the LF. The 
present sample contains 818 NB118 excess objects, 394 of which are selected as Ha emitters. Optical 
spectroscopy has been obtained for 62% of the NB118 excess objects. Empirical optical broadband 
color classification is used to sort the remainder of the sample. A comparison of the LFs constructed for 
the four individual fields covered by the observations reveals significant cosmic variance, emphasizing 
that multiple, widely separated observations are required for such analyses. The dust-corrected LF 
is well described by a Schechter function with L+ = io 43 00±0 - 52 erg s"\ = io- 3 - 20±0 - 54 Mpc~ 3 , 
and a = —1.6 ± 0.19. We compare our Ha LF and SFR density to those at z < 1, and find a rise in 
the SFR density cx (1 + z) 3A , which we attribute to significant evolution. Our Ha SFR density of 
iQ-i.oo±o.i8 Mq yr _1 Mpc -3 is consistent with UV and [O II] measurements at z ~ 1. We discuss 
how these results compare to other Ha surveys at z ~ 0.8, and find that the different methods used 
to determine survey completeness can lead to inconsistent results. This suggests that future surveys 
probing fainter luminosities are needed, and more rigorous methods of estimating the completeness 
should be adopted as standard procedure (for example, with simulations which try to simultaneously 
reproduce the observed Ha LF and equivalent width distributions). 

Subject headings: galaxies: distances and redshifts - galaxies: evolution - galaxies: luminosity func- 
tion, mass function - galaxies: photometry - galaxies: star formation 



1. INTRODUCTION 

The luminosity of the Ha nebular emission line is a 
star formation rate (SFR) indicator valued for its rela- 
tively direct physical connection to short-lived massive 
stars. In the local universe, it has been well calibrated 
(e.g., Kennicutt 1998; Kennicutt et al. 2009), and exten- 
sively used to measure the SFRs of individual galaxies, 
as well as the SFR density over cosmic volumes (e.g., 
Kennicutt & Kent 1983; Gallego et al. 1995; Salzer et al. 
2001; Gavazzi et al. 2002; Brinchmann et al. 2004; Han- 
ish et al. 2006; Meurer et al. 2006; Ly et al. 2007; Dale 
et al. 2008; Kennicutt et al. 2008; Lee et al. 2009). 

While it is desirable to extend Ha studies of galaxies 
to earlier cosmic times, such work is observationally dif- 
ficult because Ha is redshifted into the infrared beyond 
z ~ 0.4. Early attempts yielded samples of Ha emitters 
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that were small in size and did not sample representative 
cosmic volumes, because of the limited depth and areal 
coverage of the observations (e.g., Yan et al. 1999; Hop- 
kins et al. 2000; Tresse et al. 2002). Other SFR indica- 
tors, which are accessible in the optical at higher redshift 
(e.g., the rest-frame UV continuum or other bluer emis- 
sion lines), have therefore been more commonly used. 
Measurements of the SFR density have now been made 
at redshifts as high as ~ 6. A sharp, order of magni- 
tude decline is seen in the star formation activity of the 
universe from z ~ 1 to the present day (see e.g., Hop- 
kins 2004). However, the amalgam of measurements that 
constitute our current understanding of the cosmic star 
formation history shows a scatter of at least a factor of 
a few, which considerably reduces their usefulness as a 
constraint on models of galaxy evolution. 

In this context, it is important to determine the ex- 
tent to which systematics between different SFR indica- 
tors (such as the variable impact of dust attenuation and 
dependence on metallicity) contribute to the observed 
scatter in the cosmic star formation history. In partic- 
ular, it is useful to trace the history with a consistent 
indicator throughout time, and then compare the over- 
all histories determined with different indicators. This 
paper represents a step in this process, and aids in the 
robust extension of Ha measurements of the SFR density 
to higher redshift. 
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Here, we present Ha luminosity functions (LFs) and 
SFR densities at z = 0.80. At this rcdshift, the age 
of the universe is 6.6 Gyr, assuming a f2jv/, hjo] 

= [0.7, 0.3, 1.0] cosmology, which we adopt throughout. 
Our analysis is based on narrowband (NB) observations 
from the NEWFIRM Ha (NewHa) Survey (J. C. Lee 
et al. 2011, in preparation), which has been conducted 
with the NOAO Extremely Wide-Field Infrared Imager 
(NEWFIRM; Probst et al. 2004, 2008) at the KPNO 
4 m telescope. The current sample contains ~400 Ha 
emitting galaxies above 3<r significance, which have been 
identified over an area of 0.82 deg 2 . 

This paper is organized as follows. In Section 2, we 
give an overview of the NewHa data that are used in our 
analysis. In Section 3, we describe the selection of NB118 
excess emitters. We also discuss the techniques (i.e., ded- 
icated follow-up spectroscopy and empirical broadband 
color classification) that are used to identify the emis- 
sion line(s) responsible for the narrowband excess. A 
description of how we compute emission-line fluxes, lu- 
minosities, and equivalent widths (EWs) is provided in 
Section 4. In Section 5, estimates of the survey's com- 
pleteness and the surveyed volume are presented, and 
Section 6 presents the Ha LF and the comoving SFR 
density at z ~ 0.8. Comparisons of these results with ex- 
isting Ha measurements are described in Section 7. We 
also compare our results with that of other recent z ~ 0.8 
narrowband Ha surveys, and examine the reasons why 
inconsistent results may arise between these studies. A 
discussion of our results and their implications for the 
evolution of typical galaxies at z ~ 0.8 are provided in 
Section 8, and concluding remarks are provided in Sec- 
tion 9. Magnitudes are reported on the AB system (Okc 
1974). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. The NewHa Survey 

The NewHa Survey is an ongoing campaign designed 
to extend deep, wide searches for emission-line galaxies 
into the intermediate redshift universe. The survey takes 
advantage of the 27.'6 x 27'.6 field of view of the NEW- 
FIRM, which achieved first light on the KPNO 4m tele- 
scope in 2007 February, and became available for general 
observing in 2007 November. 

NewHa uses narrowband observations to identify 
emission-line galaxy candidates. Imaging is taken 
through 1% filters that arc designed to sample low OH 
airglow windows in the near-infrared, and objects are se- 
lected by detection of a photometric excess in a narrow- 
band, relative to continuum measurements in a broad- 
band. Bandpasses are centered at 1.18 /jm and 2.09 fmx, 
which capture Ha emission at redshifts of 0.80 and 2.19, 
respectively. The continuum flux is constrained with J 
and K s imaging. A combination of techniques, including 
dedicated follow-up spectroscopy and empirical broad- 
band color classification, is used to identify the emis- 
sion line(s) responsible for the narrowband excess, as 
discussed in more detail below. A full description of the 
overall survey is provided in J. C. Lee ct al. (2011, in 
preparation). Here, we give a summary of the 1.18 /im 
imaging observations, data reduction, and selection tech- 
nique used to construct a robust sample of Ha emitters 
at z = 0.80. 



2.2. NEWFIRM Observations 

The analysis presented in this paper is based upon 
NEWFIRM J band 9 and 1.18 fim narrowband (here- 
after NB118) 10 observations of areas in the Subaru-XMM 
Deep Survey (SXDS; Furusawa et al. 2008) and Cosmic 
Evolution Survey (COSMOS; Scoville et al. 2007) extra- 
galactic deep fields. The observations were carried out 
in 2007 December, 2008 September, and 2008 October at 
the KPNO 4m telescope. Data were acquired for three 
regions in the SXDS and one region in COSMOS, each 
spanning the ~ 750 arcmin 2 subtended by the detector 
array of the camera. In total, these observations cover 
0.82 deg 2 for a comoving volume of 9.12 x 10 4 /if 3 Mpc 3 
at z = 0.8. 

We follow standard near-infrared observing procedures 
to obtain the data. Integration times of 240 and 30 s 
are used for the individual NB118 and J exposures, re- 
spectively. The NB118 exposures are read-out with eight 
Fowler samples. In 2008, on board co-adding was enabled 
in the camera, and was used for our J-band observations, 
with every two J exposures being summed. A combina- 
tion of 9-, 6-, and 4-point dither patterns is used to cover 
a 75" square grid with 61 positions. The dithers serve 
to smooth over cosmetic defects, bridge the 35" gaps be- 
tween NEWFIRM's four 2048 x 2048 InSb arrays, and 
enable the rejection of pixels affected by persistent after- 
images caused by the latency properties of the detectors. 
The pattern, with small offsets from the initial position, 
is repeated as necessary to achieve a minimum 3tr depth 
of 23.5 AB (NB118) and 23.7 AB (J), in apertures con- 
taining at least ~80% of the flux of a point source, given 
the seeing conditions during the observations (see Sec- 
tion 3.1). Cumulative integration times range between 
8.2 and 12.7 hr in NB118 and 2.3 and 4.0 hr in J. The 
median seeing during our observations was ~ 1"2, and 
varied between l'.'O and 1"9, hence point sources are ad- 
equately sampled with 0"4 pixels. Sky conditions were 
mostly photometric, but some data were taken through 
thin cirrus. A summary of the observations is given in 
Table 1. 

2.3. NEWFIRM Data Reduction 

Data reduction is performed with our own automated 
PyRAF-bascd pipeline, which builds upon routines from 
the IRAF/nfextern package and is optimized for the 
processing of NewHa observations. Our procedures fol- 
low standard, iterative, near-infrared reduction tech- 
niques to produce flat-fields, subtract the sky back- 
ground, and reject artifacts, particularly those due to 
persistent afterimages from bright sources. 

The NewHa pipeline processes the data in two passes. 
One of the purposes of the initial pass is to create a deep 
object mask that is used to construct flat-fields and sky 
frames from the science images themselves. The pipeline 
first subtracts the dark current, masks for bad pixels, 
corrects for the nonlinear response of the detector, and 
performs a preliminary sky-subtraction using a median 
of the temporally closest five exposures. The geometric 
distortion is rectified, and the astrometry is calibrated 
relative to the Two Micron All Sky Survey (2MASS) 

9 A = 1.250 (im; 8\ = 0.180 /Ltm 

10 A = 1.184 Aim; <5A = 0.011 Jim 
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TABLE 1 
Summary of NEWFIRM Imaging 



Field 


R A V\nr 

rv,.rt.., ucc 


Observation Dates 


Filter 


Int. Time 


FWHM 


Limitin 


g Magnitude (AB, 3cr) 




(J2000) 






(hr) 


(arcsec) 


SE 


SW 


NE 


NW 


COSMOS 


10:01, +02:01 


2007 Dec 4,5 


J 


2.30 


1.20 


23.63 


23.71 


23.74 


23.74 






2007 Dec 2-5 


NB118 


8.16 


1.50 


23.48 


23.54 


23.55 


23.59 


SXDS-N 


02:18, -04:38 


2008 Sep 28,29, Oct 1,22,23 


J 


3.52 


1.10 


23.91 


24.03 


24.03 


24.05 






2008 Sep 29, Oct 1,22,23 


NB118 


8.47 


1.20 


23.77 


23.94 


23.90 


23.98 


SXDS-S 


02:18, -05:15 


2007 Dec 3-5 


J 


2.40 


1.25 


23.67 


23.75 


23.81 


23.78 






2007 Dec 2-5 


NB118 


10.28 


1.60 


23.40 


23.51 


23.58 


23.61 


SXDS-W 


02:16, -04:57 


2008 Sep 23-26,28, Oct 21,22 


J 


3.97 


1.20 


24.02 


24.14 


24.14 


24.19 






2008 Sep 23-26,28, Oct 21,22 


NB118 


12.67 


1.15 


23.88 


23.96 


24.06 


24.03 


Note. — 


Limiting magnitud 


es are reported for apertures sizes of 2'.'0 diameter for the SXDS-N and SXDS-W, 


and 2" 5 


for the 



SXDS-S and COSMOS. These apertures are chosen to contain a minimum of 80% of the flux of a point source for each pair of J and 
NB118 imaging, given the seeing conditions during the observations. 



catalog. The dithered science images are then projected 
onto a common pixel grid and stacked. A deep object 
mask is made from this initial stack, and applied to the 
individual science exposures on the original pixel grid. 

The first pass sky-subtracted frames are also used to 
identify artifacts due to persistent afterimages. Object 
masks are created for each individual frame, and pixels in 
a given image are compared to corresponding pixels from 
previous frames. Pixels containing object flux in consec- 
utive frames are masked. Non-science images, such as 
short exposures used to check and adjust the telescope 
pointing, are also included in this process. This method 
flags roughly 95% of visually apparent persistent arti- 
facts. The remaining 5% are faint and decay quickly 
enough in subsequent frames that they are not detectable 
in the final mosaics after the stacking of the entire data 
set. 

Flat-fields for each night are made by combining J- 
band science images, which have been masked for both 
real objects and persistent sources. The flats are used 
to normalize the response within each detector in both 
the J and NB118 images. Corrections that account for 
sensitivity variation between detectors are also applied. 

In the second pass, all the data are flat-fielded, and 
the sky subtraction is again performed, but with stacks 
of temporally neighboring frames which have now been 
masked of all sources. Any problematic frames (e.g., 
due to read-out problems or jumps in the telescope 
pointing) are rejected. The flat-fielded, sky-subtracted, 
persistence-masked frames are then combined to produce 
the final mosaic. The astrometry of the mosaics is tested 
against 2MASS sources, as well as against sources in the 
COSMOS optical catalog, and is found to be accurate to 
within 0'.'15-0'.'20 with negligible systematic offsets. 

Photometric calibration of the mosaics is performed us- 
ing 150-300 unsaturated 2MASS (Skrutskie et al. 2006) 
sources in each field. We check for systematic errors as 
a function of both J/NB118 magnitude and radius from 
the mosaic center, and find no significant offsets. The re- 
sultant zcropoints are accurate to within 0.05 mag. Ab- 
solute flux calibration for the 2MASS catalog is based on 
Vega (Kurucz 1979), hence we convert to AB magnitudes 
by convolving the Alter bandpasses with the Vega spec- 
trum, and find that to(AB) - m(Vega) = 0.87 (NB118) 
and 0.95 (J). 

3. SAMPLE SELECTION 
3.1. Selection of Narrowband Excess Emitters 



Sources that show a significant J— NB118 color excess 
are selected as emission-line object candidates. Our se- 
lection procedure follows general techniques commonly 
used in narrowband surveys (e.g., Fujita et al. 2003; Ly 
et al. 2007; Shioya et al. 2008; Villar et al. 2008; Sobral 
et al. 2009). 

To construct our sample, we first use SExtractor (ver- 
sion 2.5.0; Bertin & Arnouts 1996) in dual-image mode to 
generate source catalogs for each field. That is, sources 
are identified on the NB118 image, and fluxes are mea- 
sured on both J and NB118 images in matched apertures 
at the position of every 3cr NB118 detection. Detection in 
the J band is not required for inclusion in the candidate 
list. 

In each field, photometry is performed in two sets of 
apertures. The sizes of the apertures for each pair of 
J and NB118 images are chosen to contain at least 80% 
and 99% of light from a point source, given the size of the 
point-spread function (Table 1) and assuming a Gaussian 
profile. In the SXDS-N and SXDS-W fields, the seeing 
during our observations was between l'.'l and l'/2, and 
thus 2" and 3" diameter apertures are used. The aper- 
ture sizes are increased to 2'.'5 and 4" for the SXDS-S 
and COSMOS fields, where the seeing was worse and 
varied from 1'.'25 to l'/6. One of the purposes of using 
two aperture sizes is to allow the selection to be sensi- 
tive to galaxies with concentrated star formation as well 
as those with more extended nebular emission. In order 
to derive global quantities, however, the Ha luminosi- 
ties and EWs for all selected sources are calculated from 
fluxes measured in the larger aperture. 

Sources are considered to have a significant narrow- 
band excess if: 

1. They have a color above a minimum threshold given 
by 

A(J-NB118) > 0.2 mag, (1) 

where A(J — NB118) is defined below, and 0.2 mag 
roughly corresponds to the 5cr scatter in A( J — NB118) 
for bright point sources with 18 < NB118 < 20 mag. It is 
imposed to help exclude bright foreground sources with 
blue continua from the candidate list. 

2. The color is significant at the 3cr level 



A(J-NB118) > 3Jaj + a 



'NB118 



(2) 



Other recent near-infrared narrowband surveys have 
adopted a lower threshold of 2.5cr (Villar et al. 2008; 
Sobral et al. 2009, hereafter V08 and S09). We 
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Fig. 1. — Left: standard broadband-narrowband color magnitude selection diagrams, which illustrate the dependence of the mean value 
of J — NB118 for pure continuum sources "(J — NB118) C " on the slope of the continuum, as measured by z' — J. The points are color coded 
to correspond to different ranges of z' — J, as indicated in the plot. Sources with redder continua have lower (J — NB118) C compared with 
those with bluer continua. Right: the color magnitude diagram after the dependence on z' — J has been removed, and (J — NB118) C is 
subtracted. The resultant A( J — NB118) is the color (or narrowband) excess used to compute emission-line fluxes and equivalent widths. 
(A color version of this figure is available in the online journal.) 



also extract a secondary sample with this relaxed sig- 
nificance criterion to facilitate comparisons with these 
studies (Section 7). 

Two issues critical to the robust application of these 
criteria for the selection of line emitters are the determi- 
nation of A( J — NB118), and the calculation of accurate 
photometric errors. We discuss each of these in turn. 

Flux excess in the NB118 filter cannot be directly cal- 
culated from the J - NB118 color since the NB118 filter 
does not fall at the center of the J bandpass, but rather 
was designed to sample the low-OH airglow window to- 
ward the blue edge. This not only results in a non-zero 
mean color for pure continuum sources (J— NB118) C , but 
also leads to variation about the mean color as a func- 
tion of continuum slope. This systematic is illustrated in 
the left panel of Figure 1, which shows a standard color- 
magnitude selection diagram for one of our fields. The 
points in the diagram are color-coded to correspond to 
different ranges in the z' — J color, which provides a mea- 
sure of the continuum slope. Sources with red continua 
have lower values of (J — NB118) C compared with those 
with blue continua, as expected. Here, the J photom- 
etry is from our NEWFIRM observations, while the z' 
measurements are from Subaru/Suprimc-Cam (see Fu- 
rusawa et al. 2008 for SXDS and Taniguchi et al. 2007 
for COSMOS). 

To account for this systematic, we define 

A(J - NB118) = (J - NB118) - (J - NB118) C (3) 
where 

(J — NB118) C = f(z' — J). (4) 

A linear fit is performed for sources with J — NB118 
within ±0.25 mag of the overall mean value. We find that 
(J - NB118) C = -0.15(z' - J) + const, for z' - J < 0.5. 
For z' — J > 0.5, there is no apparent slope so a constant 
value is assumed for redder sources. Accordingly, we 
compute corrections and apply them relative to the value 
of ( J — NB118) C for sources with an intermediate z' — J of 



0.25 mag. The right panel of Figure 1 shows the color- 
magnitude diagram after this correction is applied and 
demonstrates the removal of the systematic. 

Photometric uncertainties are calculated by combining 
the errors generated by SExtractor with measurements 
of the background noise taken through a large number 
of apertures randomly placed on the images. This hy- 
brid scheme ensures that our uncertainties account for 
(1) correlated (non-Poissonian) noise arising from pixel 
interpolation during astrometric reprojection in the im- 
age reduction (captured in the random aperture measure- 
ments) and (2) local variations in the noise, for example, 
due to non-uniformity in the exposure map (captured by 
the SExtractor errors). A detailed discussion of these 
issues is given in Gawiser et al. (2006). 

The global, average uncertainty due to the background 
can be robustly measured using a large number of aper- 
tures (identical in size to the apertures used for the pho- 
tometry) that randomly sample the sky in the images. 
We perform this exercise separately for different quad- 
rants in the image to account for quantum efficiency dif- 
ferences between the NEWFIRM detectors. 11 The distri- 
bution of measurements is well fitted by a Gaussian func- 
tion, and the standard deviation given by the fit provides 
the la uncertainty. The average depths of our observa- 
tions determined in this way are given in Table 1. 

To incorporate information on local variations in the 
noise that is included in the SExtractor errors, the SEx- 
tractor estimates arc compared with those computed 
from the random apertures. We find that SExtractor 
yields lower uncertainties, which is expected since the 
program assumes that the noise is purely Poissonian and 
simply computes the error from the global pixel-to-pixel 
rms. To account for the contribution of non-Poissonian 

11 In Table 1, note that while three detectors have comparable 
sensitivities, the fourth is less sensitive by about 0.3 mag. The 
primary cause is an extra layer of anti-reflection coating, which 
was inadvertently applied to the detector. 
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Fig. 2. — Color-magnitude diagrams illustrating the selection of NB118 excess emitters in the SXDS-S, where the photometry has been 
measured in the larger choice of aperture (4"). Note that the color excess A(J — NB118) (see Equations (l)-(2)) rather than raw J — NB118 
color is plotted. The horizontal lines indicate ±0.2 mag from zero excess; the minimum accepted A(J — NB118) is 0.2 mag. The vertical 
lines indicate the 3ct depth of the NB118 photometry. The three curves shown correspond to the average values of the color excess at 3.0a, 
2.5ct, and 2. Oct significance, as determined from measurement of the image background with random apertures, and Gaussian fitting of the 
resultant distribution. Sources with A(J — NBff8) values significant at the 3ct level are highlighted with red circles, while those significant 
at the 2.5ct level are enclosed in black diamonds. The selection is done on a per-quadrant basis, to account for the varying sensitivities of 
the detectors. 



noise components, we scale the median SExtractor error 
(as a function of the object flux and on a per-quadrant 
basis) to match the average uncertainty determined from 
the random apertures. The errors from SExtractor are 
increased by 2%-51%, depending upon the aperture size 
used. These scaled SExtractor errors are adopted to eval- 
uate the significance of the color excess detected in our 
observations. 

The selection criteria are illustrated in Figure 2, which 
shows the color-magnitude diagrams for one of our fields, 
where the photometry has been measured in the larger 
choice of aperture. Horizontal lines indicating the mini- 
mum accepted color excess arc plotted, along with curves 
representing the average values of the color excess at 
3. Oct, 2.5ct, and 2. Oct significance (as determined using 
random apertures). Objects selected at 3ct (red cir- 
cles) and 2.5ct (black diamonds) are both shown. Here- 
after, quantities related to the 2.5ct samples are quoted 
in brackets unless otherwise noted. 



The selection is performed separately for the photome- 
try determined through the two sets of apertures, as well 
as for the individual quadrants in each field. The use 
of a second, larger detection aperture results in a ^10% 
increase of the number of narrowband excess candidates. 
In total, we find 150-300 [250-450] candidates per field. 
All candidates are inspected by visual examination of 
the NEWFIRM J and NB118 data, alongside publically 
available deep z' imaging. The inspection process leads 
to the removal of four sources from the sample (two are 
artifacts and two are significantly blended with adjacent 
candidates) . The overall process yields a total sample of 
818 [1218] NB118 excess emitters in the combined survey 
area of 0.82 deg 2 . A summary of the number of sources 
identified in each field is given in Table 2. A catalog 
of NB118 excess emitters, including emission-line fluxes 
and EWs, will be provided in a forthcoming paper (J. C. 
Lee et al. 2011, in preparation). 
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TABLE 2 

Summary of NB118 Selected Samples 



Field 


Area 


WNB1187 


/spec 


N Ha 


/spec (Ha) 


COSMOS 

SXDS-N 

SXDS-S 

SXDS-W 

Totals 


732.5 
745.2 
741.7 
746.6 
2966.0 


157 [253] 
157 [255] 
201 [265] 
303 [445] 
818[1218] 


65% [48%] 
59% [42%] 
75% [64%] 
52% [38%] 
62% [46%] 


59 [ 90] 
63 [ 93] 
130 [154] 
142 [185] 
394 [522] 


46% [32%] 
63% [49%] 
76% [70%] 
61% [48%] 
64% [52%] 



Note. — Surveyed area listed in units of arcmin . ./Vnbii87 gives the 
number of sources selected as emission-line galaxy candidates which have 3(7 
[2.5c] A(J-NB118) excess and are above A(J-NB118) = 0.2 mag. / 3pcc 
indicates the percentage of Nnbii87 with follow-up optical spectroscopy. 
N H a gives the number of emission-line galaxy candidates identified as Hct 
emitters by empirical color selection or spectroscopy, as discussed in Section 
3.2. /spec (Ha) indicates the percentage of Nna that has been spectroscop- 
ically confirmed. 



3.2. Identification of Ha Emitters 

An excess in the narrowband can arise from various 
emission lines that are redshiftcd into the bandpass of 
the filter. The main emission lines detected in our NB118 
observations are as follows: 

S III] AA9052,9532 at z = 0.31 and 0.24, 
S II] AA6717,6731 at z = 0.76, 
Ha at z = 0.80, 

[O III] AA4959,5007 at z = 1.39 and 1.36, 
H/3 at z = 1.44, and 
[O II] A3727 at z = 2.2. 
To isolate the Ha emitters in our sample, we use a com- 
bination of techniques. For 62% [46%] of the 3ct [2.5a] 
selected objects, spectroscopic data are available from 
our own targeted follow-up of NewHa narrowband excess 
sources (see below) with a small contribution from public 
spectroscopic datasets. For the remainder of the sample, 
a classification scheme based on optical broadband col- 
ors is used. The classification is empirically calibrated 
with a combination of (1) the subset of NB118 excess 
emitters for which spectroscopy and optical broadband 
photometry are available, and (2) galaxies in the COS- 
MOS spectroscopic catalog whose redshifts would cause 
an emission line to appear in the NB118 bandpass. 

Overall, we find that 48% of our 3cr narrowband excess 
sample (394/818) can be identified as Ha emitters. The 
fraction is slightly lower (43% or 522/1218) for the 2.5cr 
sample. Among our four fields, the fraction ranges from 
35% to 65%. This observed field-to-field variation is dis- 
cussed in the context of predictions for cosmic variance 
in Section 6.1. A summary of these statistics is given in 
Table 2. 

3.2.1. Spectroscopy of NewHct Narrowband Excess Sources 

We carried out spectroscopy of NewHa narrowband ex- 
cess sources between 2008-2009 with the Inamori Mag- 
ellan Areal Camera and Spectrograph (IMACS; Dressier 
et al. 2006) at the 6.5m Magellan I telescope. IMACS 
enables multi-object spectroscopy with slit masks over 
a 2 7. '4 diameter area (an excellent match to the field of 
view of NEWFIRM), and has good sensitivity to ^9500 
A. These two characteristics make IMACS an ideal in- 
strument for optical spectroscopic follow-up of NewHa 
NB118 excess sources, and in particular, Ha emitters at 
z = 0.80. Our observational setup yields spectral cov- 
erage from 6300 A to 9600 A with 9 A resolution, and 
captures the complete ensemble of strong rest- frame opti- 



cal emission lines blueward of Ha at z = 0.80, from [O II] 
A3737 at 6720 A to [O III] A5007 at 9030 A. The spectra 
have a median integration time of ~3.5 hr on average, 
reaching 5cr sensitivities of ~ 2 x 10~ 17 erg s _1 cm -2 . 
IMACS spectroscopy has been obtained for 56% [42%] of 
our candidates, and the data are currently being used to 
probe the attenuation and mctallicity properties of the 
sample (I. Momcheva et al. 2011, in preparation). A 
more complete description of the NewHa IMACS follow- 
up campaign is provided in J. Lee et al. (2011, in prepa- 
ration) and I. Momcheva et al. (2011, in preparation). 

Additional spectroscopy is available from observations 
carried out by the SXDS (M. Akiyama, 2008, private 
communication) and zCOSMOS (DR2; Lilly et al. 2007) 
survey teams. The sample sizes of the spectroscopic cat- 
alogs are 11975 for COSMOS and 4231 for SXDS, and 
among these, 2690 and 2386 fall in our survey regions. 
The zCOSMOS survey uses the VIMOS spectrograph on 
the Very Large Telescope 8 m, and targets I < 22.5 mag 
sources, as well as galaxies that are color-selected to have 
1.4 < z < 3.0. The existing SXDS spectroscopic catalog 
is formed from the composite of data from a number of in- 
dividual observing programs targeting a range of SXDS 
sub-samples, using various instruments, and achieving 
different depths. The mean and let dispersion of the 
magnitude distribution for the available SXDS z spcc cat- 
alog are i?c = 21.7 and 2.0 mag, respectively. Both of 
the resultant SXDS and zCOSMOS spectroscopic cata- 
logs arc dominated by objects that arc bright relative to 
our NB118 excess emitters, and hence the overlap with 
our sample is small (5% or 58/1218). 

In total, 504 [566] sources in our NB118 excess sample 
have spectroscopic redshifts. The redshift distribution 
for these galaxies is shown in Figure 3. There are 253 
[272] galaxies with redshifts between 0.790 and 0.817, 
which places Ha in the NB118 bandpass, and thus con- 
firms them as Ha emitters. 

3.2.2. Empirical Optical Color Classification 

For the remaining 38% [54%] of the NB118 excess 
sources that do not have spectroscopic redshifts, a clas- 
sification scheme based on publically available Suprimc- 
Cam broadband photometry is developed to separate Ha 
emitters from other sources in the sample. Our classifi- 
cation uses Rc i' z' in the SXDS fields, and r' i' z' in the 
COSMOS field. These sets of bandpasses arc chosen be- 
cause at z = 0.8, the Balmer and 4000 A breaks occur 
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Fig. 3. — Redshift distributions for sources selected as NB118 
excess emitters as determined from optical spectroscopy (black his- 
togram) and from photometric redshifts (gray shaded histogram). 
The redshifts of the primary emission lines detected in our NB118 
observations are labeled. (A color version of this figure is available 
in the online journal.) 

at -6500 A (in the R band), so NB118 Ha emitters will 
appear much redder in R—i' for a given i' — z' than other 
narrowband excess sources (Figure 4). 

In the SXDS fields, we examine the distribution of 
NB118 excess sources that have spectroscopy from our 
follow-up IMACS observations (Figure 4, left panel), and 
identify Ha emitters as those sources with 



R c - i' > 0.25, 
Rc-i'> 0.846(i' - z') + 0.204, 



iff 
if %' - 



- z' < 0.35 (5) 
z 1 > 0.35. (6) 



In the COSMOS field, our IMACS spectroscopic 
dataset sparsely samples the color-color diagram. Thus, 
to improve coverage of the diagram, we also include 
galaxies in the zCOSMOS spectroscopic catalog which 
lie outside the area of our imaging, but whose redshifts 
would cause an emission line to appear in the NB118 
bandpass (open triangles in Figure 4, right panel). Sim- 
ilar color criteria are adopted: 

Rc-i' > 0.39, if i'-z' < 0.30 (7) 
R c -i' > 1.08(i' - z 1 ) + 0.162, if i' - z' > 0.30. (8) 

An additional 141 [250] NB118 excess sources lacking 
spectroscopic redshifts are identified as Ha candidates 
with these criteria in both fields. 

Of course, this color method is a blunt selection tool 
compared with the use of spectroscopic redshifts — it 
leads to the inclusion of some interlopers, while missing 
true Ha emitters that do not lie within the color selec- 
tion region, as evident from Figure 4. In particular, the 
color criteria cannot distinguish between Ha and [S II] 
emitters because the wavelengths of the emission lines 
are separated by —100 A. 

Fortunately, we can estimate the contamination and 
miss rates by examining the color classification of the 
NB118 excess sources with spectroscopic redshifts. In 



the SXDS fields, 15% of the excess sources with available 
spectroscopy that lie in the Ha selection region are in- 
correctly classified as Ha emitters, and 5% are confirmed 
Ha emitters that fall outside the color selection region. 
The corresponding numbers for the COSMOS field arc 
29% and 8%. Applying these rates to the number of Ha 
emitters that are added to the Ha sample via the color 
method, we find that 25 [45] sources are potential inter- 
lopers and 8 [15] Ha emitters could be missed. However, 
the overall contamination and miss rates are far more 
limited with a high level of spectroscopic completeness 
of the sample: the estimated number of interlopers and 
missed Ha emitters are only 6% [2%] and 9% [3%] of the 
total Ha sample. 

4. CALCULATION OF Ha LUMINOSITIES 

4.1. Observed Line Fluxes, EWs and Luminosities 

Emission-line fluxes and EWs for the z = 0.8 Ha 
galaxy sample are calculated as follows. Flux densi- 
ties (erg s _1 cm -2 A -1 ) may be written as /nb = 
fc + F L /ANB and fj = fc + F L /AJ, where f c is the 
continuum flux density, Fl is the emission-line flux (erg 
s _1 cm -2 ), and ANB and AJ are the full width at half 
maximum of the NB118 (110 A) and J (1786 A) filters. 
Therefore, the emission-line flux and the continuum flux 
density are computed as 



F L = ANB 



fc 



1-(ANB/AJ) ' 

/j-/wb(ANB/AJ) 
1-(ANB/AJ) 



and 



(9) 
(10) 



By definition, the observed EW, which is a factor of 
1 + z larger than the rest-frame EW, is the ratio of the 
emission-line flux to the continuum flux density, and thus 



EW obs = ^ = ANB — 
Jc J J 



/nb - fj 



/nb(ANB/AJ)' (U) 

With this equation, the minimum A(J - NB118) = 0.2 
mag excess corresponds to an observed EW of 40 A. 
To compute luminosities, we adopt a distance of 5041 
hjQ Mpc, which corresponds to z = 0.803, the median 
redshift of our spectroscopically confirmed Ha emitters. 
Distributions of the observed emission-line fluxes and 
EWs for our sample of Ha emitters are shown in Fig- 
ure 5. 

4.2. Derived Ha Luminosities 

To determine the intrinsic Ha luminosity from the ob- 
served line luminosity, corrections must be applied for 
contamination of the flux by the adjacent [N II] lines, and 
for attenuation by dust internal to the galaxy. The cor- 
rections that we adopt arc based on standard, empirically 
calibrated relationships, which describe average correc- 
tions for ensemble populations. They should yield lumi- 
nosities appropriate for computing averaged/integrated 
quantities such as the SFR volume density. However, it 
should be noted that the corrections will not be accurate 
for individual galaxies, since the scatter in the relation- 
ships is large, as described below. 

[N II] contamination. The NB118 bandpass is wide 
enough to include flux from the [N II] AA6548,6583 12 



12 Hereafter, 
lines. 



; '[N II]" refers to both nitrogen nebular emission 
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Fig. 4. — Left: Rc —i' and i' — z' colors of NB118 excess emitters in the three fields observed in the SXDS (black filled circles). Open 
diamonds arc over-plotted on objects with spectroscopic redshifts, and are color-coded to indicate the emission line responsible for the 
narrowband excess, as described in the legend. Right: analogous plot for the one field observed in COSMOS, but with r' instead of Rc 
photometry. To improve coverage of the diagram in this field, we also include galaxies from the zCOSMOS catalog whose redshifts would 
cause an emission line to appear in the NB118 bandpass (open triangles), but lie outside the area of our imaging. The black lines in both 
panels illustrate our adopted color criteria for the selection of Ha emitters in the NB118 excess sample (see Equations (5)-(8)). (A color 
version of this figure is available in the online journal.) 
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Fig. 5. — Distribution of the observed Ha + [N II] emission-line 
fluxes (left) and EWs (right) for all identified Ha emitters, com- 
puted from photometry measured in the larger choice of apertures, 
as discussed in Section 3.1. The EWs are a factor of 1 + z larger 
than the rest-frame EWs. 

emission lines for narrowband excess selected Ha emit- 
ters. 

In the local universe, integrated spectroscopic surveys 
have found that the Ha/[N II] flux ratio is 2.3 for typ- 
ical galaxies (Kennicutt 1992; Gallego et al. 1997), 
and past narrowband surveys have often adopted a fixed 
correction of 2.3 for all selected Ha emitters. The flux 
ratio, however, has been shown to increase with larger 
emission- line EW (see, e.g., V08) and with decreasing B- 
band luminosities (Kennicutt et al. 2008; Lee et al. 2009). 
Such correlations are likely a consequence of the mass- 
metallicity relationship (e.g., Lee et al. 2004; Tremonti 
et al. 2004) . In a recent deep optical narrowband survey 
with Subaru Suprime-Cam, Ly et al. (2007) adopted 4.66 
based on optical spectroscopic follow-up of their emit- 
ters. The sample of Ly et al. (2007) consists of galaxies 
with fainter luminosities, hence are more metal-poor and 
would have a higher Ha/[N II] flux ratio on average. 

For this study, we follow V08 and S09 in adopting an 



EW-dcpcndent Ha/[N II] flux ratio to facilitate compar- 
isons of results. The EW-dependcnt correction was con- 
structed from thousands of z ~ 0.1 star- forming galaxies 
from the Sloan Digital Sky Survey fourth data release 
with which V08 determined the mean relationship be- 
tween the rest-frame EW of Ha + [N II] A6583 and the 
Ha/[N II] flux ratio. The correlation exhibits a scat- 
ter of ~0.2 dex, which implies that estimates of the 
Ha/[N II] ratio for individual sources will only be accu- 
rate to ~50%. Such errors will average out in the calcu- 
lation of integrated quantities however, and are adequate 
for our purposes here. 

Also, our correction assumes that the Ha/[N II] rela- 
tion does not evolve with redshift. This is a valid assump- 
tion since the evolution of mctallicity between z = 0.07 
and z = 0.7 has been found to be no more than 0.1 dex 
(see, e.g., Mannucci et al. 2009, and references therein). 
Near-infrared multi-object spectroscopy will be needed 
to examine the Ha/[N II] ratio for NB selected galaxies 
on a case-by-case basis. Assuming the V08 [N II] correc- 
tion, the Ha/[N II] ratio varies from 1.85 to 10.0 with a 
median (average) of 2.50 (2.81) for our population of Ha 
emitting galaxies. 

Dust attenuation. To correct for dust attenuation, we 
adopt a luminosity-dependent extinction relation follow- 
ing Hopkins et al. (2001): 



log[SFR obs (Ha)] = log[SFR : 
x log 



0.797 log [SFRj 
2A 



(Ha)] - 2.360 

(Ha)]+3.786 



(12) 



Hopkins et al. (2001) showed that attenuations computed 
from this equation show a 0.2 dex scatter relative to those 
based on Balmer decrement measurements for individual 
galaxies. Much like the [N II] correction, these dust ex- 
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tinction corrections are more reliable when considering 
multiple sources of a given luminosity /SFR. The Hopkins 
et al. (2001) relation was also adopted in other studies 
(e.g., Ly et al. 2007; V08; Dale et al. 2010). For our 
sample, this correction (A[Ha]) ranges from 0.66 to 1.87 
mag with a median and average of 1.19 and 1.21 mag, 
respectively. 

One piece of evidence that supports the Hopkins et al. 
(2001) relation is the work of Garn et al. (2010), which 
has looked at the UV, Ha, and mid- infrared fluxes for 
a sample of narrow-band selected Ha emitters to de- 
termine dust attenuation. They find a correlation be- 
tween the observed Ha luminosity and dust attenuation 
that is similar to Hopkins et al. (2001). This multi- 
wavelength comparison will be conducted for our sample 
of Ha emitters. We will also investigate dust attenuation 
determined from Balmer decrements, which are obtained 
from the combination of our NB118 measurements and 
our IMACS follow-up spectroscopy (I. Momcheva et al. 
2011, in preparation). Preliminary results do indicate 
that the attenuations based on Balmer decrements for 
the NB118-selected Ha emitters are consistent with the 
Hopkins et al. (2001) correction. 

5. COMPLETENESS OF THE SURVEY 

5.1. Monte Carlo Simulations 

Motivation. To account for the detection limits and 
photometric selection of the NcwHa survey, we need to 
estimate the completeness fraction as a function of lumi- 
nosity, k(L). The greatest uncertainty for such a task is 
the intrinsic EW distribution, that is, the inherent distri- 
bution of the population that we are able to seek. Since 
the narrowband and broadband imaging sensitivities of 
the survey are known, the completeness is determined for 
a source of a given brightness with a J — NB118 excess 
color (i.e., an EW). However, sources of a given emission- 
line luminosity can be faint with large EWs or bright 
with low EWs. Thus, the adopted EW distribution af- 
fects the estimated completeness for individual emission- 
line luminosity bins. For this reason, we conduct Monte 
Carlo realizations of our data, and follow a "maximum 
likelihood" approach to determine the completeness of 
NcwHa. 

Technique/approach. The methodology intends to si- 
multaneously reproduce the shape of the observed 13 
Ha+[N II] EW distribution and the observed Ha LF. We 
begin with a grid of models that adopt certain EW distri- 
butions. With these distributions, we generate artificial 
emission-line galaxies. We add noise to these galaxies, 
and use their measured magnitudes and colors to deter- 
mine if these mock galaxies satisfy the NB118 excess se- 
lection. We construct the observed EW distribution and 
Ha LF from the sample of mock NB118 excess emitters 
to compare to the respective distributions from NewHa. 
We determine a best model and define the completeness, 
k(L), as the ratio of the number of mock NB118 excess 
emitters to the total number of generated artificial galax- 
ies: 

k{L) = N(L) obs /N(L) tot . (13) 

13 Throughout this section, wc use "observed" to denote that the 
EW distributions arc a factor of 1 + z larger than the rest-frame 
EW distributions. 



Assumptions. For the prior distributions, we make two 
assumptions. First, for the Ha LF, a proper normaliza- 
tion of bright-to-faint sources is necessary, otherwise we 
will generate relatively too many bright sources, and a 
decline in the simulated LF will not been seen and can- 
not be compare to the observed Ha LF. To constrain the 
ratio of bright-to- faint sources, we examine the galaxy 
number counts (N[J]) in the NewHa fields and find that 
log [N( J)] oc 0.344 x J. This relation implies that the 
number of galaxies at magnitude J' + 1 is a factor of 2.2 
more numerous than those at J'. 

Second, each plausible model assumes a log-normal 
EW distribution described by a mean (log(EWo/A)) 
and a sigma a [log (EW/ A)]. We adopt such a distri- 
bution, since it is fairly similar to the observed NewHa 
distributions and the distributions seen locally (Lee et 
al. 2007). Our models span (log(EW /A)) = 1.15- 
1.55 and <r[log (EW/A)] = 0.15-0.40 both in increments 
of 0.05 dex for a total of 54 EW distributions. Here, 
(log(EWo/A)) is the logarithm of the rest-frame EW, 
and we later include a factor of 1+z for the "observed" 
magnitudes, colors, fluxes, and EWs. 
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Fig. 6.— Color-magnitude diagram (A(J-NB118) vs. NB118) 
illustrating the selection of mock NB118 Ha excess emitters in the 
SXDS-W for the NE detector. We show modeled galaxies as black 
filled squares, and the gray points are the repeated measurements 
of these galaxies with photometric noise included. This figure rep- 
resents 2,004,800 mock galaxies with every fiftieth point shown 
for lower resolution. The horizontal line is the minimum 0.2 mag 
excess that we adopt, and the three black curves show the la, 
2.5cr, and 3cr A(J — NB118) color excesses. The model assumes 
(log(EW /A)) = 1.15 and er [log (EW/A)] = 0.15. 

Implementation. For each of the 54 models, we gen- 
erate ~40,000 mock galaxies with J-band (continuum) 
magnitudes chosen to follow the above relation between 
log [N( J)] and J. We assume that these galaxies are spa- 
tially unresolved and follow the FWHM of seeing for 
each mosaic. For the EWs, wc randomly draw from 
a log-normal distribution given by (log(EW /A)} and 
er[log (EW/A)]. With these two information, we "fill" the 
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color-magnitude diagram (A(J — NB118) versus NB118; 
see Figure 6). 

We add observational uncertainties to the J and 
NB118 magnitudes as follows. For each mock galaxy, 
we have a randomly chosen position in the mosaics, and 
generate 200 realizations based on the sensitivity at that 
position. The sensitivity is governed by the exposure 
maps for both the NB118 and J-band mosaics, thus ac- 
counting for detector and pixel-to-pixel variations. We 
compute the NB118 or J-band 3a limiting magnitudes, 
miim{x,y), signal-to-noise ratios, S/N(x, y), and la un- 
certainties, Am(x,y), as 

mum(x,y)= (mum) + 2.5 log (^/t(x, y)/toj , (14) 
S/N(a;,j/) = 3 x io-°- 4 [™-™H m (^)] j an d (15) 



Am{x, y) 



2.5 log 



1 



S/N(x,y) 



(16) 



where (mum) and to are the individual detector's median 
limiting magnitudes and exposure times (see Table 1), 
respectively, and m is the NB118 or J magnitude. For 
the A(J — NB118) colors, the errors are added in quadra- 
ture: \J AraNBiis^, y) 2 + Amj(x, y) 2 . The magnitudes 
for the 200 realizations follow a Gaussian distribution 
with a mean of m and la of Am(x, y). 

We identify mock galaxies that are NB118 excess emit- 
ters using the same methods described in Section 3.1. 
For the smaller apertures (2" or 2'.'5), we assume that 
the flux enclosed is between 81% and 86% of the total 
flux within the larger apertures (3" or 4"). We deter- 
mine these values by considering a Gaussian distribution 
with a FWHM equal to the seeing size for each NEW- 
FIRM pointing. Note that we adopt limiting magnitudes 
appropriate for the size of the measurement aperture. 

Comparisons with observations. For each model, we 
construct the observed EW distribution and Ha LF from 
the mock NB118 excess emitters. Since we generate 
many mock galaxies, we normalize each distribution to 
best match the observed NewHa distributions. We il- 
lustrate an example of these model-observation compar- 
isons in Figure 7 for the SXDS-W field. We find the 
best fit by minimizing % 2 through a comparison of the 
mock and the NewHa EW distributions and Ha LFs. 
We show in Figure 7 that the best-fit model is able to 
reproduce both the EW distribution and the Ha LF. We 
find that observed EW distribution provides a better con- 
straint on the best-fit model (versus the Ha LF). That 
is, multiple EW models can sufficiently explain the Ha 
LF. We find the model with (log(EW /A)) = 1.35 and 

<r[log (EW/A)] = 0.40 best matches the observed distri- 
butions. 

We calculate the completeness as the ratio of the num- 
ber of mock galaxies that meets the NB118 excess selec- 
tion criteria to all mock galaxies (N w 200 x 40, 000) as 
a function of L(Ha). We illustrate in Figure 8 how the 
completeness differs by adopting different (log (EWo/A)) 
values. This comparison shows that we miss bright galax- 
ies with low EWs, which is not surprising given the 
minimum observed EW of 40 A. This fi gurc also re- 
veals that the completeness can vary between 60% and 
90% for observed fluxes of 5 x 10 _17 -6 x 10 -16 erg s _1 
cm~ 2 or roughly 0.1-1 J* (see Section 6.1), and empha- 
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Fig. 7. — Comparisons between the modeled predictions and ob- 
servations of the EW distribution (top) and the Ho galaxy num- 
ber counts (bottom) for the SXDS-W field. The shaded and un- 
shaded histograms are for the SXDS-W and all four NewHa fields, 
respectively. The dashed line in the top panel is the minimum 
adopted J — NB118 excess of 0.2 mag or an observed EW of 40 
A. The model (black curves) assumes (log(EWo/A)) = 1.35 and 
erpog (EW/A)] = 0.40, and is scaled to match the observed distri- 
butions. This model is able to produce the shape of the observed 
EW distribution and the Hct number counts. 



sizes that proper choice in the prior EW distribution is 
necessary for an accurate determination of the LF. We 
also show the completeness as a function of Ha lumi- 
nosity and observed EW for each of the fields assuming 
the best-fitting model with (log(EW /A)} = 1.35 and 
ct [log (EW/A)] = 0.40 in Figure 9. 

Discussion of simplifying assumptions. In the above 
calculations, we made a couple of simplifications to al- 
low the Monte Carlo simulations to be significantly less 
computationally intensive. First, we assume that the 
galaxies are unresolved sources. Second, instead of de- 
tecting sources on the images, we estimate the measured 
magnitudes and fluxes where the photometric uncertain- 
ties are based on pixel-to-pixel and detector-dependent 
sensitivity. Third, it is important to note that our cur- 
rent simulations assume that the EW distribution is the 
same at all continuum magnitudes. It would be more 
correct to allow the EW distribution to vary as a func- 
tion of continuum magnitude, since it has been shown 
that more luminous star-forming galaxies tend to have 
lower emission-line EWs than those at lower luminosity 
(e.g., Lee et al. 2007). However, this adds another de- 
gree of freedom to our models, and significantly increases 
the computational time required to complete the simu- 
lations. Future work will include greater complexity in 
the modeling of the EW distribution, and will examine 
the impact of the current assumptions on the resultant 
completeness corrections. These simplifications allow for 
a few orders of magnitude faster production of complete- 
ness results rather than the approach of adding sources 
directly to the images, and is able to reproduce the distri- 
butions of scatter in Figures 1 and 2 (e.g., see Figure 6). 
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Fig. 8. — Completeness correction derived from Monte Carlo sim- 
ulations of our SXDS-W data for models with crpog (EW/A)] = 
0.40. The Ha fluxes have not been corrected for any dust extinc- 
tion. The color and line style conventions are for different adopted 
(log (EWo/A)), as indicated in the lower right. This figure illus- 
trates that as a larger (log (EWo/A)) is adopted, the complete- 
ness increases for bright emission-line galaxies. As expected, the 
adopted minimum J — NB118 color excess of 0.2 mag is unable 
to identify galaxies that are very bright with low EWs. (A color 
version of this figure is available in the online journal.) 

Prior to these MC simulations (for this discussion, we 
refer to the above approach as "Method II"), we gener- 
ated a simulation where artificial extended sources were 
directly added to the mosaics and detected. In this sim- 
ulation (referred to as "Method I"), we adopted a log- 
normal Gaussian EW distribution with (log (EWo/A)) = 
1.52 and c[log (EW/A)] = 0.16, motivated by the results 
in Lee et al. (2007). We performed Method I by adding 
1000 artificial galaxies to each NEWFIRM mosaic, and 
compared the number of artificial NB118 excess emitting 
galaxies to the number detected. We created the galaxies 
using IRAF/mkobjects with physical parameters similar 
to those found in local galaxies. The parameters include 
luminosities, Ha EWs, semimajor to semiminor axial ra- 
tios (between 0.15 and 1.0), and the Ha disk scale length 
(3.6 kpc; Dale et al. 1999). After extracting sources in 
the same manner described in Section 3, we calculated 
the completeness. We repeated each simulation 30 times 
using different seed numbers to increase the statistical 
accuracy of the results (N ~ 30, 000 per NEWFIRM 
pointing) . The similarities in the shape of the complete- 
ness curves in Methods I and II suggest that the simpli- 
fications we made do not have a significant effect on the 
estimated completeness. There are certain differences 
such that the comparison is not completely "apples-to- 
apples" ; however, good agreement suggests that proceed- 
ing with Method II is satisfactory. The greatest benefit 
of this simulation is the determination of a model that 
best matches two different sets of measurements. 

Finally, we assume a power-law distribution for N(J) 
to normalize the number of bright and faint sources. Ide- 
ally, the distribution should be constructed from galaxies 



at z ~ 0.8. An investigation with the COSMOS photo- z 
sample (after the simulation was completed) finds that 
the faint-end slope of the number counts is approximately 
0.25, which is to be compared to the adopted 0.34. Also, 
an exponential decline exists at bright continuum lumi- 
nosities. This will impact our completeness estimates 
for the brightest emitters. Future work will address and 
examine the impact that these differences have on the 
completeness reported here. 

5.2. Effective Volume 

The volume that the survey is capable of probing is 
determined by the shape and width of the NB filter. The 
ideal filter will have a perfect top-hat profile, and will 
survey the same co-moving volume at all emission-line 
fluxes. However, with real filters a weak emission line 
can either result from the line falling in the wings of 
the NB filter profile or from an intrinsically weak line 
located near the filter center. The net result is that a 
weak emission line is more likely to be detected near filter 
center, so a non-square filter reduces the effective volume 
at the faint end (Ly et al. 2007). 

We can quantify how much the filter profile deviates 
from being a perfect top-hat by comparing the area un- 
derneath the profile with that of a rectangle with width 
equal to the FWHM, and height equal to the maximum 
transmission of the actual filter. The area enclosed by 
the profile is just 9% smaller, so the shape of the band- 
pass is close to the ideal. 

To determine the effective volume as a function of 
emission-line flux, we calculate the range in wavelengths 
such that an emission line is considered detectable within 
the NB filter. This emission line has an intrinsic S/N. We 
then place it at different wavelengths to determine what 
the degradation in the S/N will be due to lower through- 
put. We define an emission line to be undetected below 
2.5<7, and this criterion yields the minimum and maxi- 
mum NB118 wavelengths (hence redshift) that is observ- 
able for a particular S/N. The effective comoving volume 
per unit steradian would then be 
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Here, Z\ and Z2 refer to the minimum and maximum 
rcdshifts that the emission line is detectable. Note that 
these equations assume a flat universe. The maximum 
surveyed volume (V c s/Q = 1.11 x 10 5 Mpc 3 deg~ 2 or 
A A = 110 A) is observable for S/N > 5.1 and decreases 
to V cS /n = 5.54 x 10 4 Mpc 3 dog" 2 (AA = 55 A) at 
S/N = 2.57. This is illustrated in Figure 10. 

6. RESULTS 

6.1. Luminosity Function and Star Formation Rate 
Densities 

Figure 11 presents the (observed) Ha LF for this sur- 
vey with and without completeness corrections. Apply- 
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Fig. 9. — Completeness corrections as a function of luminosity (k(L); left) and EW (k(EW); right) derived from Monte Carlo simulations 
of our data. The line types denote different fields: COSMOS (solid), SXDS-N (dotted), SXDS-S (dashed), and SXDS-W (dot-dashed). The 
best model ((log (EWo/A)) = 1.35 and cr[log (EW/A)] = 0.40) is used. The higher completeness for SXDS-W and SXDS-N is due to higher 
spatial resolution and sensitivity. The vertical dashed line in the right panel is for the minimum observed EW limit of 40 A. For k(EW), 



we consider mock galaxies with an Ha emission-line flux above the 2.5a flux limit of sb 1.2 X 10 
erg s" 1 cm -2 when the adjacent [N II] emission is included (i.e., the total NB118 emission-line flux). 



erg s * cm 



" 2 which is » 1.8 X lO -17 




The median variation of the number density relative 
to the average of all four fields is ^50%, and we illus- 
trate in Figure 13 the fluctuation of the four different 
pointings relative to the average. Using predictions from 
Somerville et al. (2004, hereafter S04), we find that the 
expected fluctuation per NEWFIRM pointing (shaded 
region in Figure 13) is consistent with what is observed. 
The hourglass-like shape of expected amount of field-to- 
ficld fluctuations is due to (1) the stronger clustering of 
luminous galaxies, (2) the decrease of cosmic variance 
with number density, and (3) the small volume surveyed 
at low luminosities due to the weakness of emission lines. 
Thus, the minimum of ficld-to-field variations is around 
an extinction-corrected Ha luminosity of 1 x 10 42 erg s" 1 . 

It is common to model the Ha LF by fitting the LF 
with the Schechter (1976) function: 

, / L \ a ( L\dL 

$>(L)dL 



L 



— ) CX P(-— ) T - (21) 



S/N Ratio 

Fig. 10. — Effective surveyed volume. The FWHM (in wave- 
length) of the surveyed volume versus S/N ratio. The horizon- 
tal line corresponds to 50% of the maximum effective volume 
(V eS /fl = 1.11 X 10 5 Mpc 3 deg~ 2 ). 

ing the necessary completeness corrections and adopting 
the Hopkins et al. (2001) extinction correction, we have 
the Ha LF shown in Figure 12. Both raw, and extinction- 
and completeness-corrected number densities as a func- 
tion of luminosity are listed in Table 3. The binned LF 
can be summed together to obtain a model-free, lower- 



limit Ha luminosity density of C = 10 



40.01±0.04 



MpC" 3 [C : 



^q40.08±0.03 er g g -l 



Mpc" 



erg 



This function was derived from the Press & Schechter 
(1974) formalism, which describes the halo mass func- 
tion, and was adapted to explain the distribution of 
galaxy continuum luminosities. An explanation for why 
the Schechter function can describe the continuum LF 
is that a connection exists between the luminosities of 
galaxies and their stellar masses and/or halo masses. 
However, such a link may be much weaker with the Ha 
luminosity and SFR; thus, one might question whether 
the Schechter function is a good model for the Ha LF. For 
the current analysis we simply assume that the Schechter 
function can be used to adequately model the distribu- 
tion of Ha luminosities and SFR since it appears to pro- 
vide a good fit to our data. This also facilitates com- 
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TABLE 3 
Ha Luminosity Function atz~ 0.81 









2.5ct 








3a 




logL 


N 


-/Vspcc 


*(L) 






-/Vspcc 


4>(L) 










Raw number densities 








4U. /U 


11 





11.17±3.4 


0.08 










a n on 


60 


8 


49.34±6.4 


0.22 


20 


3 


1 O 7/1 _LO Q 


U.zZ 


a a in 
41. 1U 


116 


27 


83.76±7.8 


0.47 


65 


19 


4U.oO±D.U 


U.4y 


A A QH 

41. oU 


102 


50 


63.02±6.2 


0.72 


78 


45 


4o.yt)±c>.z 


n 7Q 


41. OU 


98 


74 


55.39±5.6 


0.79 


96 


73 


O4.Ul±D.0 


n 70 


41. /U 


66 


51 


36.21±4.5 


0.79 


66 


51 


OO.Zl±4.0 


u. <y 


/i i on 
41. yU 


34 


26 


18.65±3.2 


0.81 


34 


26 


lO.ODzto.Z 


n q 1 


/t o 1 n 
4z. 1U 


19 


18 


10.42±2.4 


0.80 


19 


18 


in ^ nil) a 
1U.4Z±Z.4 


n en 
U.oU 


/to Qn 


8 


7 


4.39±1.6 


0.81 


8 


7 


a Qn_Li a 


U.ol 


/i o ten 


8 


8 


4.39±1.6 


0.86 


8 


8 


A Q O J_ i a 


U.ou 




Extinctior 


i and completeness-corrected number 


densities 




40.90 


3 





59.79±34.5 


0.07 










41.10 


18 


1 


245.40±57.8 


0.14 


3 





43.48±25.1 


0.18 


41.30 


55 


8 


156.58±21.1 


0.28 


19 


1 


46.46±10.7 


0.26 


41.50 


77 


15 


108.82±12.4 


0.52 


13 


10 


50.12±7.6 


0.53 


41.70 


99 


11 


93.94±9.4 


0.69 


64 


37 


54.58±6.8 


0.71 


41.90 


89 


55 


66.39±7.0 


0.78 


84 


53 


61.75±6.7 


0.78 


42.10 


63 


18 


43.66±5.5 


0.80 


63 


18 


43.66±5.5 


0.80 


42.30 


50 


13 


34.85±4.9 


0.79 


50 


43 


34.85±4.9 


0.79 


42.50 


29 


21 


19.70±3.7 


0.81 


29 


21 


19.70±3.7 


0.81 


42.70 


18 


14 


12.32±2.9 


0.80 


18 


14 


12.32±2.9 


0.80 


42.90 


12 


12 


8.12±2.3 


0.81 


12 


12 


8.12±2.3 


0.81 


43.10 


5 


1 


3.37±1.5 


0.81 


5 


4 


3.37±1.5 


0.81 


43.30 


4 


4 


2.46±1.2 


0.89 


4 


1 


2.46±1.2 


0.89 



Note. — is normalized to 1 x 10 4 Mpc 3 dcx 1 and luminosities (L) are given 

in erg s — . k(L) is the survey completeness fraction defined in Section 5.1. Numbers 
reported in the top half arc prior to any completeness corrections while completeness is 
included for the bottom half. The spectroscopic completeness (A^ spcc ) as a function of 
luminosity is shown. 



TABLE 4 

SCHECHTER FITS, Ha LUMINOSITY DENSITIES, AND SFR DENSITIES 



Survey 


logL* 


log 
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log£(0) 


log£(L lim ) 


log Psfr(0) 


log PSFR(0) a 


log PSFR(ilim) a 


NewHa; 
NewHa 
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S09 
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42.33«; 1 J 
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40.01 ± 0.08 
40.04 ± 0.05 
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-1.00 ± 0.18 
-0.98 ±0.05 
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-0.96 


-1.10 ±0.08 
-1.06 ±0.05 



Note. — and <&* in units of erg s - and Mpc~ , respectively. Luminosity (£) and SFR densities (psfr) are provided for L > and 
L > Lii m . All LFs adopt Hopkins et al. (2001) dust attenuation. Note that S09 incorrectly reported the normalization for their Hct LF 

and was off by a factor of ~0.43. D. Sobral provided a new preliminary fit, which adopts the Hopkins et al. (2001) extinction and the proper 
LF normalization. 
1 



Corrections for AGN contamination applied. 

parisons with previous work. However, this assumption 
should be explored further when more accurate LFs from 
future studies show evidence that a Schechter function is 
not the best model to explain the Ha LF. 

In order to obtain the best-fitting Schechter parame- 
ters, a Monte Carlo simulation was performed to con- 
sider the full range of scatter in the extinction- and 
completeness-corrected Ha LF. Wc ignored luminosities 
below our 50% completeness limit (1.0 x 10 41 erg s _1 
observed; 3.0 x 10 41 erg s _1 extinction-corrected). Each 
datapoint was perturbed randomly 1 x 10 times follow- 
ing a Gaussian distribution with ler in given by 
Poisson statistics. Each iteration is then fitted to obtain 
the Schechter parameters. The best-fitting Schechter 
parameters are then determined from the averages of 
these fits. The confidence contours for the best fit 
are shown in Figure 14. A summary of the best fits 
and the corresponding integrated Ha luminosity density 



(£ = J L^(L)dL) and SFR density (see below) is pro- 
vided in Table 4. We find a relatively steep faint-end 
slope for the Ha LF (a = —1.6) at z ~ 0.81, indi- 
cating that galaxies below 0.2 and 1L± contribute 61% 
and 89% to the total Ha luminosity/ SFR density, re- 
spectively. Often, past studies have opted to fixed the 
faint-end slope since they were unable to reliably con- 
strain it. Following this methodology, wc also report the 
results of our fits when a is set to —1.6 in Table 4. 

The extinction-corrected Ha luminosity density can be 
converted into a SFR density by using the recipe given in 
Kennicutt (1998): SFR(Ha) = 7.9 x 10" 42 i(Ha) where 
the SFR is given in M@ yr _1 and the Ha luminosity is 
given in erg s . This conversion assumes a Salpctcr IMF 
with minimum and maximum masses of 0.1 Mq and 100 
Mq and solar metallicity. We determined that the Ha 
SFR density is p SFR = io-°- 96 ± 018±0 04 M Q yr" 1 Mpc" 3 
down to L = 0, where the second sets of errors account 
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Fig. 11. — Observed (i.e., no dust attenuation corrections have been applied) Ha LF before (left) and after (right) applying completeness 
correction. Corrections have also been made for [N II] contamination. The combined measurements from the four NewHa fields are shown 
as filled squares, where as open squares represent the individual pointings. V08 and S09 measurements are shown as asterisks and triangles, 
respectively. The axis range is kept the same as other plots of the LF for comparison purposes. All LFs have selected sources down to 
2.5(7 significance. Conversion between Ha luminosities and SFRs is based on the Kennicutt (1998) relation. (A color version of this figure 
is available in the online journal.) 



for cosmic variance estimated from S04. Compared to 
measurements at z < 0.1 (Gallego et al. 1995; Perez- 
Gonzalez et al. 2003; Brinchmann et al. 2004; Nakamura 
et al. 2004; Hanish et al. 2006; Ly et al. 2007; Westra et 
al. 2010), our Ha SFR density at z ~ 0.8 is higher by a 
factor of 3.8 to 16.6 with a median of 8.1. 

It is thought that the Ha luminosity density is dom- 
inated by emission from star formation and not active 
galactic nuclei (AGNs). However, to accurately compute 
the SFR volume density, we statistically account for the 
fraction of the Ha luminosity volume density originat- 
ing from AGNs. While estimating the AGN fraction 
is an observational challenge because deep x-ray data 
and/or spectroscopic information are needed, previous 
studies have typically found that 10%-15% of galaxies 
have AGNs. For example, Brinchmann et al. (2004) es- 
timate that 11% of the Ha luminosity density is due to 
AGNs at z < 0.1. Gallego et al. (1995) found 15% AGN 
contamination to the Ha SFR density at z ~ 0. V08 
used X-ray data for a small (~50) sample of Ha emit- 
ters at z ~ 0.8 and found 10%±3%. S09 looked at the 
[O III]/H/3 and [O II]/H/3 flux ratios (Rola et al. 1997), 
for a subset of 28 z ~ 0.8 Ha emitters with optical spec- 
troscopy and reported 15%±8%. Our preliminary analy- 
sis of the rest-frame optical emission-line flux ratios from 
our IMACS spectroscopy (see Section 3.2.1) finds simi- 
lar results. Based on the [O III]/H/3 and [O II]/H/3 flux 
ratios, 5 (10) of 141 Ha emitters are AGNs (LINERs). 
Thus, we correct the above Ha SFR density by 11% to ac- 
count for AGN and LINER contamination. This reduces 
our total Ha SFR density to p SFR = i -i-°o±o.i8±o.04 
M Q yr" 1 Mpc -3 . 

7. COMPARISONS WITH OTHER NEAR-INFRARED Ha 
STUDIES 



Recently, two other independent groups have per- 
formed relatively wide-field near-infrared narrowband 
imaging on 3-4 m class telescopes to search for high- 
z emission-line galaxies. V08 first surveyed 626 arcmin 2 
for Ha emitting galaxies at z ~ 0.84. They identified 165 
galaxies and obtained a extinction-corrected SFR den- 
sity of 0.17 M Q yr" 1 Mpc" 3 . S09 surveyed a total of 1.3 
deg 2 , identified 743 Ha emitting galaxies at z ~ 0.84, 
and determined a SFR density of 0.1 M Q yr -1 Mpc -3 . 
NewHa complements these surveys through a combina- 
tion of depth and volume surveyed: it covers almost 5 
times more area than V08, and while S09 covers about 
50% more area, our survey is 0.6 dex deeper. These 
advantages simultaneously allow us to (1) obtain better 
constraints on the faint-end slope and "knee" of the LF 
(we acquired ~10% accuracy on the slope and 0.5 dex 
on L+\ see above), and to (2) reduce ficld-to-field fluc- 
tuations. The NewHa dataset not only enables us to 
compute a more robust estimate of the LF, it also allows 
us to better understand the properties of sub-L* galaxies 
and their role in the overall star formation history of the 
universe. Comparisons of the LFs between these three 
surveys are provided in Figures 11, 12, and 14 and Ta- 
ble 4. We begin by comparing how NB excess emitters 
arc selected and how Ha emitters arc identified, and then 
discuss the discrepancies between these surveys that are 
apparent. 

Redshift and sensitivity. The V08 and S09 surveys 
probed volumes at a slightly higher redshift of ~ 0.84. 
Of course, significant evolution between z ~ 0.81 and 
z ~ 0.84 (corresponding to At = 127 Myr at ~ 6.5 Gyr) 
is not expected. The main consequence of the different 
surveyed redshift arises from impact of sky background 
level at the wavelength of redshiftcd Ha. NewHa tar- 
gets a cleaner window in the sky spectrum. This partly 
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Fig. 12. — Extinction-corrected Ha LF at z ~ 0.81. Dust extinction corrections adopted the Hopkins et al. (2001) equation. The color 
and point-style conventions follow those in Figure 11. The left figure shows each of our NEWFIRM pointing and the average with the 
best-fitting Schechter function (L* = 10 43 00±0 - 52 erg s _1 , = l()- 3 ' 20±0 - 54 Mpc~ 3 , and a = -1.6 ±0.19) as the solid line. On the right 
panel, we compare our average LF with the LFs of V08 (shown as asterisks) and S09 (shown as triangles). S09 measurements were adjusted 
to adopt the Hopkins et al. (2001) dust extinction equation instead of A(Hct) = 1.0 mag. Note that the ordinate axes have different ranges. 
Conversion between Ho luminosities and SFRs is based on the Kennicutt (1998) relation. (A color version of this figure is available in the 
online journal.) 
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Fig. 13.— Field-to-field fluctuations in the Ho LF for the NEW- 
FIRM pointings. The rr-axis shows the extinction-corrected Ha 
luminosity while the y-axis shows the number density of sources 
normalized to the average. Color convention of points follows those 
used in Figure 11. The shaded regions represent the la variation 
predicted from the ACDM model of S04. (A color version of this 
figure is available in the online journal.) 

leads to the factor of two and four times deeper depth in 
emission- line flux sensitivity compared to V08 and S09. 

Selection of excess emitters. All three surveys iden- 
tify NB excess emitters above 2.5a significance in the 
J— NB color (though we also report results based on a 
more robust 3er-selectcd sample). However, there are dif- 



FlG. 14. — Confidence contours for the Schechter fit. 68% and 
95% level for L*, and a arc shown from a Monte Carlo simu- 
lation of the Ho LF with extinction, completeness, and [N II] con- 
tamination corrections. The faint-end slope is fixed to o = —1.6 
for the upper left panel and is free for the other three panels. For 
these other three panels, we overlay the Schcchtcr-fitting results of 
V08 and S09 as asterisks and triangles, respectively. 

ferences in the minimum J— NB color criterion and the 
aperture(s) used. For example, we required at least a 
A(J— NB118) color of 0.2 mag and used two aperture 
sizes (2"or 2.5" and 3 or 4"). S09 used a 3" diameter 
selection and required a minimum J — NB color of 0.3 
mag. V08 indicated that they selected sources in a to- 
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tal of 10 apertures out to 5 times the FWHM with a 
minimum J — NB color of 0.15 mag. We found that the 
inclusion of larger apertures for NewHa only provided 
10% more candidates that the smaller aperture failed to 
catch. Larger apertures would only allow V08 to identify 
bright and extended galaxies, which are often at lower 
redshifts. 

Identification. In all three surveys, spectroscopic red- 
shift is available for some NB excess emitters to dis- 
tinguish other emission lines from Ha. The follow-up 
spectroscopy that we have classifies 62% [46%] of our 3a 
[2.5a] NB excess emitter candidates. On the other hand, 
V08 and S09 classified the majority of their NB excess 
emitters with photometric redshift. These surveys have 
follow-up spectroscopy for 9% (138/1527; S09) and 48% 
(69/165; V08). It is difficult to compare directly their 
spectroscopic completeness against ours since (1) S09 did 
not report the size of the UDS spectroscopic sample, thus 
the 9% is a lower limit on their spectroscopic complete- 
ness, and (2) V08 only reported spectroscopy for their 
Ha emitters rather than the full NB excess emitter sam- 
ple. For the latter, NewHa has 52% spectroscopic com- 
pleteness of Ha emitters. Nevertheless, the spectroscopic 
completeness of NewHa is higher than those of V08 and 
S09. 

Similarities in the LF. The observed number density 
of Ha emitters is shown in Figure 11 (left panel). It 
illustrates that all three surveys agree to within ^50%, 
which indicates that the differences in the selection of 
Ha emitters and NB excess emitters do not significantly 
affect the observed number densities of galaxies. The 
V08 observed LF is generally higher, but it is consistent 
with our SXDS-S observations, which has a similar area 
coverage, so cosmic variance is likely the cause (see below 
for further discussion). 

In Figure 12 (right panel) we show the extinction- and 
completeness-corrected LF for all three Ha surveys. It 
is apparent that the bright end of S09's LF is consis- 
tent with NewHa's. This is to be expected since both 
surveys cover large areas and incompleteness corrections 
are less of an issue. In addition, V08 is in agreement 
with NewHa below a luminosity of ^10 42 erg s _1 . This 
is also expected since V08 reaches comparable sensitiv- 
ity to our survey. However, there are two discrepancies 
worth addressing. 

Differences in the LF. First, the extinction- and 
completeness-corrected LF of V08 is 0.2-0.3 dex higher 
at L*. However, the luminous end of the V08 LF is not 
well determined since their survey consists of three point- 
ings (two in the Groth strip, one in GOODS-N) totaling 
less than 0.2 deg 2 . Thus, the combination of Poisson 
fluctuations and cosmic variance can explain the higher 
number density at the bright end. 

The second discrepancy is that S09 find a higher num- 
ber density of faint Ha emitters, after survey complete- 
ness corrections arc applied (sec Figure 11, right panel). 
S09 claimed that the completeness correction is a fac- 
tor of 2-3 at their 2.5tr flux limit, while our Monte Carlo 
simulation indicates that 50% completeness occurs at our 
~ 6a flux limit. Furthermore, S09 find that their com- 
pleteness gradually decline while our Monte Carlo sim- 
ulation shows that above an observed Ha emission-line 
luminosity of ~ 10 415 erg s _1 , the completeness is gen- 



erally >80% and falls of rapidly with fainter luminosi- 
ties. How rapidly the completeness declines will affect 
the shape of the LF, and thus the completeness correc- 
tions can alter the determined Schechter parameters. An 
illustration of this is provided in Figure 14 where the dif- 
ferences in the Schechter parameters for these two studies 
are shown. Recall that the two LFs prior to any com- 
pleteness corrections are consistent. 

The differences between NewHa and S09 at the faint 
end imply that (1) deeper data must be obtain such that 
incompleteness determinations are less of an issue or (2) a 
more standardized rigorous procedure is needed for com- 
pleteness estimates. For example, a method that consid- 
ers a maximum likelihood approach that simultaneously 
produces the observed emission-line EW and LF (such 
as those performed in Section 5.1). 

8. DISCUSSION 
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Fig. 15. — SFR density from Ha surveys. Open squares are 
31 measurements from the literature (see Section 8 and Table 5 
for references), while NewHct measurements are shown as the 
filled circle (LF integrated to L = 0) and as a triangle (above 
the survey limit). Our Ha measurements have been systemati- 
cally reduced by 11% to account for potential AGN contamina- 
tion. The uncertainties in our "total" SFR density include an es- 
timate of the amount of cosmic variance expected for our survey 
(see text) and the uncertainties in fitting the LF with a Schechter 
profile. The dashed line is the fit adopted by Dale et al. (2010) 
for z < 2 measurements (this fit excluded our measurement): 
log (psfr[M yr" 1 ]) = -2.06 + 3.39 log (1 + z). 



In this section, we compare our Ha SFR density mea- 
surements with those published in the literature for a 
range of redshifts. We limit the comparison to other 
Ha-based measurements to avoid systematic issues with 
other SFR indicators. 

The latest compilation of Ha measurements was made 
by Dale et al. (2010). We add our measurements to this 
compilation and plot them as a function of redshift in 
Figure 15. All of the measurements plotted are summa- 
rized in Table 5. We note that the measurements re- 
ported here correct for a few mistakes found in the orig- 
inal papers (see Table 5 footnotes). The dashed line is 
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a fit determined by Dale et al. (2010) where the NewHa; 
SFR density was not included in the fitting process. It 

has the form of log ( M ^-i ) = -2.06 + 3.39 log (1 + z). 

Our SFR density measurement with the removal of 11% 
for AGN contamination is above the Dale et al. (2010) 
fit, but consistent within the uncertainties. This relation 
indicates that the Ha SFR density increases by a factor 
of ~ 10 per unit redshift at z < 1.5. 

To understand this redshift evolution, we compare in 
Figure 16 the confidence contours of the Schechter pa- 
rameters for measurements at z = 0.07, 0.39, and 0.81. 
The values reported for z = 0.39 is based on a com- 
bination of data from two complementary surveys: the 
Subaru Deep Field (SDF; Ly et al. 2007, extremely deep 
for 0.25 deg 2 ) and the Wyoming Survey for Ha (Dale et 
al. 2010, shallower sensitivity by 1.8 dex but covers ~1 
deg 2 ). The binned LFs are combined together and fit 
with a Schechter function. Likewise, z < 0.1 measure- 
ments from the SDF are combined with Gallego et al. 
(1995). We find that the characteristic luminosity sys- 
tematically increases by 0.95 dex (0.7 dex) from z ~ 0.1 
[z rj 0.4) to z ~ 0.8, while the normalization is similar 
at all three redshifts. This indicates that the increase in 
the SFR density is a result of z ~ 0.8 galaxies pro- 
ducing stars at a rate that is ~10 times that of local L+ 
(~ 10 42 erg s~ J ) galaxies. 

Infrared surveys have provided complementary con- 
straints on the evolution of the LF and SFR densi- 
ties. These measurements are sensitive to UV radia- 
tion of massive stars that are absorbed by dust and re- 
radiated at rest wavelengths of 10-100 /im. Lc Floc'h 
et al. (2005) examined the infrared luminosity of 24^m- 
selected galaxies and found that for z < 1, the LF 
evolves as log(L*jR) cx 3. 21^2(1 + z) and log (</>*,m) oc 

0.7lQ g(l + z). These results of significant evolu- 
tion and weak <!>* evolution are consistent with those of 
NcwHa, granted this can be for a different population of 
galaxies. 

In addition, we find that our Ha SFR density mea- 
surement at z k, 0.8 is consistent with z ~ 1 UV and 
[O II] SFR density measurements (Hopkins 2004, and 
references therein), though the scatter among the mea- 
surements spans a factor of ~ 2. The causes for the large 
scatter include (1) cosmic variance, (2) systematic issues 
involving the SFR indicators, and (3) the different selec- 
tion biases that affect each of the surveys. These issues 
also affect the LF and SFR density derived from infrared 
surveys. All of these issues will be addressed over the 
next few years to improve the accuracy that the cosmic 
star formation history can be determined. In general, 
cosmic variance will gradually become less of an issue 
with surveys covering at least several deg 2 . In future 
work, data from the NewHa survey, in combination with 
UV data from GALEX, mid-infrared fluxes from Spitzer, 
and [O II] fluxes from follow-up spectroscopy (see Section 
3.2.1), will directly enable us to compare the SFRs based 
on these indicators in hundreds of individual galaxies to 
address point (2). 

And finally, upcoming studies, including ones based 
upon the NewHa Survey, will investigate the differences 
in galaxy samples that result from the use of different 
selection techniques. For example, infrared surveys pri- 



marily probe the dustiest galaxies, while UV surveys 
preferentially select the bluest galaxies, and emission-line 
selected surveys are biased against low-EW galaxies. Di- 
rect comparisons of samples selected with different tech- 
niques in the same volumes will allow us to account for 
variations in the SFR volume density, which are due to 
such selection biases. As mentioned in the Introduction, 
ultimately it will also be possible to trace large fractions 
of the cosmic star formation history with a single indica- 
tor, and then compare the consistently measured histo- 
ries from multiple indicators. 
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Fig. 16. — Redshift evolution in Schechter parameters. Confi- 
dence contours for L*, <£>*, and a for H« measurements at z ~ 0.1 
(black circles), z = 0.39 (dashed line), and z = 0.81 (solid line). 
These contours are derived from a Monte Carlo realizations of the 
Ha LFs. The greatest difference is in L* with it increasing by ~ 1 
dex. 



9. CONCLUSIONS 

We have presented new measurements of the Ha LF 
and SFR volume density for galaxies at z ~ 0.8, based 
on 1.18/Um narrowband imaging from the NewHa Survey. 
With a 3cr Ha emission-line flux depth of « 1.9 x 10 7 
erg s _1 cm -2 (a luminosity of ~ 6 x 10 40 erg s _1 ) and 
an area coverage of 0.82 deg 2 , the NewHa survey allows 
for a reduction of field-to-ficld fluctuations to 10%, and 
for robust estimates of the faint-end slope (~10% accu- 
racy) and the location of the "knee" (0.5 dex accuracy) of 
the LF. We have identified 818 NB excess emitters above 
3cr, and 394 are classified as Ha emission-line galaxies 
at z « 0.80. The classification utilized a large spectro- 
scopic sample providing unambiguous determination of 
redshifts for 62% of the sample. These spectra were also 
used to calibrate the multi-color selection of the remain- 
ing Ha emitters without spectroscopic follow-up. 

We constructed the extinction- and completeness- 
corrected Ha LF. Corrections for [N II] flux contami- 
nation and the effective surveyed volume, as a function 
of line flux, were applied. The LF is well described by 
a Schechter function with L* = io 43 00±0 - 52 erg s _1 , 
= io- 3 - 20±0 - 54 Mpc" 3 , and a = -1.6 ± 0.19. When 
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TABLE 5 
Compilation of Ha SFR Densities 







Area a 


jy 


io & IPSFRJ 


Gallego et al. (1995) 


0.022±0.022 


471.4 deg 2 


176 


-1.91 ±0.04 


Tresse & Maddox (1998) 


0.20±0.10 


500 


138 


-1.61 ±0.03 


Yan et al. (1999) 


1.3±0.5 


~ 85 


33 


—0.574 ± 0.182 


oullivan et al. (2UU0J 


n i c _l_n i c 
U.15±U.15 




216 


— l.ou ± U.Ub 


Tresse et al. (2002) 


0.73±0.30 




30 


i r»c+0.07 

- L06 -o.os 


Pujita et al. (2003) 
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706 


348 


— l.OU_ 17 


Hippelein et al. (2003) 
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407 


92 


1 QQ + 0.10 

-1.83_0.13 


Perez-Gonzalez et al. (2003) 


0.025±0.025 




79 


± -°- L — 0.08 


Brinchmann et al. (2004) 


0.10±0.01 


SDSS 




-1.54 ±0.07 


Nakamura et al. (2004) 


0.06±0.06 
0.079±0.013 


SDSS 


1482 


_•, q 4 +0.106 

-1.87 ±0.03 


Hanish et al. (2006) 


0.06±0.06 


SINGG 


110 


-i.8oig:J? 


Ly et al. (2007) 


0.08±0.015 


868 


318 


-1.87±0.29 d 




0.24±0.011 


868 


259 


-2.11 ±0.24 d 




0.40±0.018 


868 


391 


-1.79 ± 0.20 d 


Geach et al. (2008) 


2.23±0.016 


0.60 deg 2 


55 


-1.00° 


Morioka et al. (2008) 


0.242±0.009 


875+SDSS 


575 


-1 456+°'?° 


^IViinvFi pf ftl f9nns'i 
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n 94-i-n noo. 
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V Uo 
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Westra &c Jones (zUUoJ 
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707 
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bnim et al. (__U09j 
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1 ^ Hno- 2 
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Dale et al. (2010) 


0.16±0.02 


4.19 deg 2 


211 


-2.002 ± 0.20 




0.24±0.02 


4.03 deg 2 


424 


-1.877 ±0.21 




0.32±0.02 


4.13 deg 2 


438 


-1.691 ± 0.23 




0.40±0.02 


1.11 deg 2 


91 


-1.660 ±0.25 


Hayes et al. (2010) 


2.19±0.014 


56 


55 


f 


Westra et al. (2010) 


0.05±0.05 


4 deg 2 


322 


-2.18 ±0.10 




0.15±0.05 


4 deg 2 


1127 


-1.92 ±0.09 




0.25±0.05 


4 deg 2 


1268 


-1.82 ±0.05 




0.34±0.04 


4 deg 2 


848 


-1.81 ±0.03 


NewHa (total) 


0.809±0.008 


0.82 deg 2 


522 


-1.00±0.18 d 


NewHa (L > L lim ) 


0.809±0.008 


0.82 deg 2 


414 


-1.10±0.09 d 



Unless otherwise indicated, areas are in arcmin . 
2 Psfr in units of Mq yr — 1 Mpc — 3 . Corrections for dust extinction have been included. 
These values integrated the LF to L — except for the last line in this table. 
^ See Ly et al. (2007) for discussion of potentially 50% contamination, which was not 
accounted. 

^ Estimates for cosmic variance are included within the uncertainties. 

^ We determined that the normalization of the LF, _ *, was reported incorrectly in these 
papers by a factor of ~ 0.43. Here, we report the correct SFR densities. 

Since Hayes et al. (2010) used the LF of Geach et al. (2008) for constraints on the 
bright end, and the normalization of the Geach et al. (2008) LF was reported incorrectly, 
the reported SFR density of Hayes et al. (2010) is thus affected and not reported here. 



the LF is integrated to L = 0, we determine a SFR den- 
sity of psfr = 10- 100±018±0 04 M Q yr" 1 Mpc" 3 . This 
SFR density is (on average) 8.1 times higher than z < 0.1 
measurements. We determined that the characteristic 
Ha luminosity is systematically higher at z ~ 0.81 by 
0.70 and 0.95 dex compared to z = 0.39 and z ~ 0.1 
estimates, respectively. The normalization of the LF at 
z ~ 0.8 is similar to what is seen for z ^ 0. This may 
imply that the cause of the redshift evolution in the Ha 
SFR density is a result of z ~ 0.8 L± galaxies producing 
stars at a rate that is «10 times that of typical galaxies 
seen locally. 

The depth and completeness of current high-z Ha sur- 
veys significantly limit the accuracy to which the (1) Ha 
SFR density, (2) the shape of the LF, and (3) its evolu- 
tion can be measured. This is underscored by the fact 
that three independent Ha surveys (including NewHa) 
show excellent agreement in the observed number densi- 



ties, but exhibit discrepancies at the factor of two level 
after completeness corrections are applied. These dis- 
crepancies lead to different conclusions on the evolution 
of star-forming galaxies: while NewHa shows redshift 
evolution in L+, V08 report evolution in both L+ and 
_>*, and S09 find evolution in The differences be- 
tween the results presented here and that of V08 can be 
attributed to cosmic variance and small number statis- 
tics at the luminous end. Fully understanding the differ- 
ences between our results and that of S09 require a more 
in-depth comparison of the completeness corrections for 
both studies. Future surveys probing fainter luminosi- 
ties, which will circumvent uncertain completeness cor- 
rections above emission-line fluxes of ~ 2 x 10~ 17 erg 
s _1 cm~ 2 , in combination with rigorous simulations for 
incompleteness that are consistently applied across stud- 
ies, will allow for convergence on the true cosmic SFR 
volume density, and better understand the factors that 
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drive the evolution. 

The OH background is the primary limitation for 
ground-based near-infrared surveys, which explains the 
historical dearth of Ha measurements at z = 0.5-0.8 and 
in the -ff-band window (z ~ 1.5). These epochs will be 
studied with the new Hubble/WFCS infrared grism. For 
example, Atek et al. (2010) has begun a survey, that will 
yield Ha SFRs at z = 0.25-1.6 as well as other emission 
lines (e.g., [O II] and [O III]) out to z ~ 4 to probe a 
significant fraction of the early universe. Both the NB 
and grism surveys complement one another through a 
combination of surveyed area, depth, and redshift. 

This study illustrates that Ha can be extended to high 
redshift, and with more sensitive detectors and wide field 
coverage in the near future, Ha measurements for thou- 
sands of galaxies at z « 1-3 will be possible to trace 
cosmic star formation history with a consistent SFR in- 
dicator over the past 11 billion years. 
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